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SECTION 1 
INTRODUCTION 


This report describes work performed under Contract NAS8-36039 on several 
existing solidification models for which computer codes and documentation were 
developed under Contract, NAS8-33573. The models describe r .ne solidification 
of alloys In which there is a time-varying zone of coexisting solid and liquid 
phases; i.e., the S/L zone. The primary purpose of the models is to calculate 
macrosegregation in a casting or ingot which results from flow of interden- 
dritic liquid in this S/L zone durino solidification. The flow, driven by 
solidification contractions and by gravity acting on density gradients in the 
interder.dr i t ic liquid, is modeled as flow through a porous medium. In Model 
1, the "steady-state model," the heat flow characteristics are those of 
steady-state solidification; i.e., the S/L zone is of constant width and it 
moves at a constant velocity relative to the mold (1). In Model 2, the 
"unsteady-state model," the width and rate of movement of the S/L zone are 
allowed to vary with time as it moves through the ingot (2). Each of these 
models exists in two versions. Models 1 and 2 are applicable to binary 
alloys; models 1M and 2M are apolicable to multicomponent alloys (3). 

Several enhancements to the models and their associated data bases are des- 
cribed here. Two new multicomponent alloys, MAR-M246(Hf) and PWA 1480 (Alloy 
454), were adoed to the data base. The development of the data for these 
alloys is described in Appendices A and B; sample results are in Section 2. A 
new graphical output selection lias been added to all four programs to allow 
display of the region in the S/L zone where gas-caused porosity may develop. 
This feature is described in Section 3 and in Appendix C. All programs devel- 
oped or modified under contracts NAS8 -33573 and NAS8-36039 were installed on a 
VAX* computer at Marshall Space Flight Center. The operating and maintenance 
procedures for this installation of the materials processing programs are in 
Section 4. 


*V AX is a trademark of Digital Equipment Corporation 
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SECTION 2 

NEW MULTICOMPONENT ALLOYS 


The original solidification models were designed for binary alloys (1,2). 
They operate from a data base containing all input data 'or each alloy, so 
that applying the models to a new binary alloy just requires an addition to 
the data base. Under a previous contract, parallel versions of the steady- 
state and unsteady-state rodels were developed for multicomponent alloys. The 
change from binary to multicomponent alloys required substantial model changes 
(3). The mul t icomponent models also operate from a data base, so that 
applying the model to the two ne-* alloys, MAR-M246(Hf ) and f’WA 1^80, only 
required that the alloy data be developed and included in the data base. The 
data were developed as shown in Appendices A and B by Dr. Poirier of the Uni- 
versity of Arizona. These data include the liquid and solid metal densities 
and a local approximation to the phase diagram. Due to the limited time 
available under this contract, data were estimated by a purely analytical pro- 
cess which used the published properties of related alloys and pure metals, 
rather than determined empirically. The multicomponent data base now contains 
three alloys; the alloy selection page of the interactive model input is 
shown below: 


1 MAR-M246 

2 MAR-M2^6(Hf) 

3 PWA 1^80 (ALLOY A5M 


ENTER ALLOY NUMBER 


Sample cases from the steady-state and the unsttidy-state models for both new 
alloys are in the following subsections. 
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ii 4 . li x i • ' ■ 

Or. FGCii . 4 

2.1 MAR-M246(Hf) 

The properties for MAR-M24o(Hf ) are shown in Fig. 2.1. Steady-state and 
unsteady-state casts have been run for the same input conditions as the MAR- 
M2A6 cases in Sections 3«^*3 and 3.^. 1 » of Ref. 3. The graphical output from 
those cases follows. 


REFERENCE TRm°era'lRE 
EUTECTIC T EM°ERA7ljRE 
SCUD DENSITY 
REFERENCE DENSITY 
SOLID EUTECTIC DENSITY 
DENSITY/TEMPERATURE DEPENDENCE 
VISCOSITY 
SURFACE TENSION 


I.38OE-03 DEG C 
’ .230E-03 DEG C 
7.820E*ee GM/CM*«3 
7.820E*00 GM/CM**3 
8.0S0E*C0 GM/CM*.»3 
-! . 2S0E-03 GM/CM**3 / DEG C 
4.000E-02 GM/'CM*SEC) 

I . 700E-03 DYNES/CM 


COMPONENT 

COMPOS I T I OK 
(WEIGHT 
PERCENT) 

EQUILIBRIUM 

PARTITION 

RATIO 

TEMPERATURE/ 

COMPOSITION 

DEPENDENCE 

DENS I "Y/ 

COMPOSITION 

DEPENDENCE 

HYDROGEN 
’NTERACT ION 
. JEFFICIENT 

CR 

9. 000E*G0 

9.000E-01 

-2.200E*00 

-I.300E-02 

3. G00E-03 

CO 

l.000f.*t)1 

! . I08E*00 

0.000E*e0 

1.000E-03 

3. I00E-03 

ft 

5. 5001 >00 

1 . I00E+00 

-3.200E*00 

-I . I00E-0I 

1 .400E-0? 

T I 

i.500L<«e 

e.seec-at 

-7. 708E+&Q 

-3.500E-02 

1 . 300C-02 

W 

1 . 000 E +01 

1 . 200E*00 

3.5 00ET00 

4.200E-02 

1 . I63E-02 

MO 

2.S00E+00 

B.S00E-01 

-1 .000E*00 

3. 300E-02 

1 . i00E-02 

7 A 

' . 500E+00 

9. 70eE-01 

-2. 500E+00 

3.500E-02 

1 . 7 0BE-02 


7.500E-01 

2. 500E-01 

-6. 100E+01 

-7.200E-0I 

6.500E-02 

HT 

1 . 500E *-00 

2. 300E-0I 

-1.040E+01 

5.200E-02 

! . I00E-02 


Figure 2.1. Properties of MAR-M2l+6( Hf ) 

2.1.1 STEADY-STATE EXAMPLE 

Graphs are shown from the steady-state model operating on MAR-M2A6(Hf ) for the 
same solidification conditions as in the MAR-M2A6 case in Ref. 3. Fig. 2.2 is 
a plot of volume fraction liquid across the S/L zone. Fig. 2.3 shows the 
i nterdendr i t i c fluid flow. Fig. 2. A shows the final local average weight per- 
cent of each element. The maximum fluid velocity in this case is only about 
one-fourth the maximum velocity in the MAR-M2A6 case. As in the case of MAR- 
IAS, the vectors are all vertical and downward, and there is no macrosegreaa- 
t ion. 
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Figure 2.3. I nterdendr i t ic Fluid F? 
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Figure 2.». Final Local Aver-., ._‘ght Percent of Each Element 



2.1.2 UNSTEADY-STATE EXAMPLE 

The output shown here is directly comparable to the output for MAR-M2^6 shown 
in Ref. 3. The first two plots (Fig, 2.? and 2.6) show the interder.dri t ic 
fluid flow in the S/L zone and contours of the weight percent of chromium in 
the final solid when the solidus has progressed to 3 cm from the chill and the 
liquldus is at 4 cm from the chill. The remaining plots (Fig. 2.7 through 
2 . 15 ) show the vertical profile of each element in the solidified casting. 
Because the solidification is vertical, there is essentially no segregation 
except for the inverse segregation near the chill. 
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Figure 2.5. I nterdendr i t i c Fluid Flow in the S/L Zone (Plot 1) 
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Figure 2.1*». Vertical Profile of Each Element in Sc’idified Casting (Plot 10) 
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Figure 2.15. Vertical Profile of Each Element in Solidified Casting (Plot II) 
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2.2 PWA U80 (ALLOY 4jjA) 

The properties for PWA 1 480 ore shown In Fig. 2.16. Steady-state and 
unsteady-state cases have been run for the same Input conditions as the MAR- 
M246 cases In Sections 3«^»3 and of Ref. 3. The graphical output from 
those cases are shown In Sections 2.2.1 aid 2.2.2. 


RE c ERF.fCE TEMPERATURE 
E J T EC > IC TEMPER4TJ3E 
SOulD density 
«E r ERENCfc DENSITY 
5CLID EUTECTIC DENSITY 
density/ temperature dependence 

VISCOSITY 
GLHTPCE TENSION 


I . 4Gtt£*03 DE3 C 
1.2C0E-O3 DEG C 
7.770E*0ii GM/CM*«3 
7. .'20E*00 GH/CH**3 
7. &60E*00 3M/CM*«3 
-1.250E-03 GM/CM**3 / DEG C 
4.000E-02 GH/CH*SEC> 
I.710E-03 DYNES/CM 


COMPONENT 


CP 

r.c 

ri. 
r I 

w 

TA 


COPOSITION 
(WEIGHT 
•PERCENT J 

1.0006*01 
5. 000E *00 
5.0<WE*00 
l.500E*00 
4.000E*00 
1.2OK40: 


EOL' I LIBRIUM 

PARTITION 

RATIO 

9.00PE-0I 
I. 100E*00 
Q.500E-01 
5.S00E-01 
1 . 2G0E ♦00 
7. 000E-01 


TEMPERATURE/ DENSITY/ 
COMPOSITION COMPOSITION 
DEPENDENCE DEPENDENCE 


-2.200E*00 

0.0006*00 

-3.200E*00 

-7.700E*00 

3.5006*00 

-2.5006*00 


-I . 3006-02 
1.000E-03 
-1.1006-01 
-3.500E-02 
4.200E-02 
3.500E-02 


HYDROGEN 

INTERACTION 

COEFP!CIEN T 

3.G00E-03 
3. 100E-03 
.400E-02 
1 .300E-02 
1.100E-C2 
1 . 100E-02 


Figure 2.16- Properties of PWA lA90 (Alloy 4 5 k) 
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2.2.1 STEADV-STATE EXAMPLE 

This example Is directly comparable to the steady-state example In Section 2.1 
and to the MAR-M246 example In Ref. 3. The same three plots a'e shown (Fig. 
2 . 1 7 v 2 . 1 8 , and 2.19). In the cate of PWA 1^80, th» velocity vectors near ’•he 
solidus are vertical and point upward indicating some expansion of the final 
material to solidify. In this case, the liquid eutectic has a density of 8.12 
gm./cc while the solid eutectic has a density of 8.05 qm„/cc. The maximum 
velocity occurs in the upward flowing material. U i th vertical, steady-state 
solidification, no macrosegregation Is predicted. 
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Figure 2.19. Final Local Average Weight Percent of Each Element 


2.2.2 UNSTEADY-STATE EXAMPLE 

The output shown here is the same plot selection as for MAR-M2*»6 shown in Ref. 
3 and for MAR-M2*»6(Hf ) shown in Section 2.12. The first two plots (Fig. 2.20 
and 2.21) show the i nterdendr i t ic fluid Mow in the S/L zone and contours of 
the weight percent of chromium in the final solid when the solidus has 
progressed to 3 cm from the chill and the liquidus is at cm from the chili. 
The remaining plots (Fig. 2.22 through 2.27) show the vertical profile of each 
element in the solidified casting. Because the solidification is vertical, 
there is essentially no segregation except for the inverse segregation near 
the chill. 
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Figure 2.27. Vertical Profile of Each Element in Solidified Casting (Plot 8) 


SECTION 3 

POROSITY PREDICTION 


A new selection wes added to the graphical output of all four models. The new 
plot shows the region In the S/L zone where hydrogen bubbles may nucleate. 
the steady-state models, this plot may be generated at the end of each case 
(Fig. 3.1). The plot can be repeated for any Input values of ambient pressure 
and initial amount of dissolved gas, because these values do not affect the 
solidification calculation. In the unsteady-state models, the ooroslty plot 
may be generated at any checkpoint during the solidification while the S/L 
zo.ie exists (Fig. 3.2). 

The model predicts pore formation at a point in the S/L zone if 

p.. > p ♦ p + 2 n/r (3.1) 

o 

where 

p^ - partial pressure of dissolved gas, 

p Q - ambient pressure, 

p ■ pressure due to fluid flow and meta I lostat i c height, 

0 - surface tension, and 

r - bubble radius. 

Discussions of the Inequality (3*1) and of porosity in Al-Cu alloys can be 

found in Ref. A and 5* 

Before each plot is generated, the user is asked to input the ambient pressure 
and the initial concentration of hydrogen in th*.: liquid metal (Fig. 3-1 and 
3.2). The latter is used together with data on the solubilltv of hydrogen 
in the alloy to calculate the partial pressure of dissolved hydrogen in the 
interdendri tic liquid at the given temperature and composition. The solubil- 
ity data and the liquid metal surface tension are stored in the multicom- 
ponent-alloy data base for each alloy. In the binary-alloy models, there is 

no data bese change because only aluminum-copper admits significant amounts 
of dissolved hydrogen. Calculation details of the partial pressure appear in 
Appendix C. 


3-1 



JRIGINAL PAUi- 
OF POOR QUALITY 


PWA 

1 483 -ALLOY 454) STEADY-STATE MODEL 


26- APR - 


function to plot 


PLOT TYPE 

I 

riNAu _OCAl AVERAGE COMPOSITION 

1 

Y - profiles 

2 

MACROS* GREGATI ON 

2 

X - profiles 

3 

VELOCITY 

3 

VECTOR FIELD 

4 

PRESSURE: P-P0 



5 

PRESSURE - BULK YDROSTATIC P 



6 

FRACTION LIQUID 



7 

MASS FLOW 



B 

SOLUTE FLOW 



g 

POROSiTY 




t N 7 -.R ITEM NUMBER OF FUNCTION TO PLO'. OR 
P TQ PROCEED. 


Lst.mctc of a-'factad by gas-causad porosity 

(Valle fo*- \drogan only. ) 

Enter asbiant presev.-* (atmospneresi: 1. 

Entur anoun* of dlaaol.'ad gas initially present In woiten r.atal iwt. pet.!: 5.E-4 

Figure 3*1. Selection of Porosity-Prediction Graphics in 
Steady-State Models 


' J wr. Mea (ALLOT 4i4. JNS’EADy-STATE MODEL 


15-MAv-OS 15:34:45 


' UNC T 1 CN rj PLO" 


P.-w'' T"PE 


■«•••»«•• S/L /ONE 

* VELOCITY 

? PP£S5U-»£: P-P0 

J PRESSURE - BULK •-(''DPOSTATIC P 

4 FRACTION LIQUID 

5 .IASS F_OW 

6 SOLUTE FLOW 
/ POROSITY 


1 Y PROFILES 

2 X PROFILES 

3 VECTOR FIELD 

4 CON T 0JPS 


“»*•*»< INGOT *«*•** 

8 F:r*w_ local average cobosition 
) vrjLUK fraction EUTECTIC 


FNTLR I TEN NUMBER OF FUNCTION TO PLOT. OR 
P TO PROCEED. 


! st urate o* -aglcns v'(ecc.ed by gas-caused porosity 
i»c!id for Syoroyen only, j 

Ertor ambient pnn3.*-e ' atmosp^e-es <: 1. 

E-ter amount ot disso/ad gas initially present in moltan natal (wt. pet. i: l.E-S 

Figure 3 . 2 . Selection of Fo. os i ty-Pred ict ion Graphics 
in Unsteady- State Models 


3-2 


3.1 B I NARY ALLOYS 

The only binary a’loy system for whicn the models make porosity prediction is 
aluminum-copper. Appendix E shows the development of the equations for the 
partial pressure of hydrogen dissolved in the i nterdendr i t ic liquiu of solidi- 
fying aluminum-copper alloys. Appendix F presents the calculation of the sur- 
face tension of the Interdendritic liquid. The constant coefficients from 
these two calculations were not put into the binary-alloy data base because 
the form of the equations probably would not generalize to ofher bi.,ary alloy 
systems which exhibit gas-caused porosity. Aluminum-copper is the only alley 
system currently in the binary-alloy data base which admits significant 
amounts of dissolved hydrogen is aluminum-copper. 

Figures 3.3 through 3-5 show the predicted regions of hydrogen bubble nuclee- 
cion in A1-4.5 wt.pct. Cu undergoing horizontal steady-state solidification. 
They also illustrate the effect of the ambient pressure and the initial con- 
centration of dissolved gas on the extent of the predicted region. In Fig. 
3*3» gas bubble nucleation is predicted over about 40 percent of the S/L zone 
when the ambient pressure is 1 atm and the initial concentration of hydrogen 
is 2 x 10"^ wt. pet. In Mg. 3.4, all condition® are the same except that the 
ambient pressure has been reduced to 10 ^ atm. The model now predicts bubble 
nucleation over about 90 percent of the S/L zone. A similar effect is seen in 
Pig. 3*5 where the ambient pressure is 1 atm, but the initial concentration of 
hydrogen has been increased tc 2 x 10 ^ wt. pet. 

3 . 2 MULT1C0MP0NFNT ALLOYS 

The form of the partial pressure calculation is the same for al! multicom- 
ponent alloys currently in the data base. The surface tension and hyd.ooen 
interaction data have been added to tne data base for each alloy. Porosity 
predictions for a case of unsteady-state sol idi f ; cat ion of PWA 148C follow. 
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3.2.1 MODIFICATION TO DATA BASE 

The surface tension, o , and the hydrogen interaction coefficient for each 
element, e™ have been added to the data base. The structure of the 
multicomponent alloy data base is now: 



Card 3 is repeated for each constituent in the alloy. Cards 1 through 3 are 
repeated for each alloy. 

3.2.2 PWA 1480 EXAMPLE 

Figures 3-6 through 3-10 show predictions of regions of hydrogen bubble 
nucleation in the time-varying S/L zone of PWA 1480 undergoing vertical 
directional unsteady-state solidification subject to the same input conditions 
as in the example in Section 2.2.2. In Fig. 3*6, the S/L zone has just begun 

to form at the chill surface, the liquidus isotherm has only progressed to 
0.02 cm from the chill, and the fraction solid at the chill is still small. 
For an ambient pressure of 1 atm and an initial hydrogen concentration of 

8 x 10 wt. pet., the model predicts no bubble nucleation in the S/L zone at 
this point in the solidification. Later, after the liquidus isotherm has 

progressed to 0.4 cm from the chill, but with some liquid still remaining next 
to the chill, the mode 1 predicts bubble nucleation over about 40 percent of 
the S/L region (Fig. 3 • 7 ) - Figures 3*8 through 3*10 show the porosity predic- 
tions when the liquidus isotherm is at 0.8, 2.0, and 4.0 cm from the chill. 

There is no variation in the Y direction, because for directional solidifica- 
tion solidification, the velocity vectors are all vertical; hence, the 

pressure field is iniform in the Y direction. 
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SECTION 4 

VAX- INSTALLATION OPERATING PROCEDURES 
AND MAINTENANCE GUIDE 


4.1 VAX/VMS SYSTEM USAGE 

The materials processing programs have been installed on a Digital Equipment 
Corporation (DEC) VAX* 11/780 computer at Marshall Space Flight Center. The 
VAX/VMS operating system is an extremely user-friendly environment; the models 
can be maintained with little specialized knowledge of the system. This 
document is intended to help users run the models and make mino- modifications 
of the code. It is recommended that users gain a more comprehensive knowledge 
of the VAX/VMS system from the available DEC documentation before undertaking 
extensive pregram modifications or program development. 


4.1.1 SOME COMMANDS 

A few rudimentary commands are described in Table 4.1. A good introduction to 
elementary system use is in the Primer listed in Section 1.2. 


Table 4.1. Rudimentary Commands 


Command 

Descr ipt ion 

DELETE f i lename. typ;n 

Deletes version n of file f i lename. typ. 

DIRECTORY 

Lists the files in the current default directory. 
Use the HELP command to learn the DIRECTORY com- 
mand qualifiers that can provide more information 
such as file size or creation date. 

EDIT f i lename. typ 

Invokes the system text editor. After you are in 
the editor, you can get help from a HELP command. 
Use th? (EOT) editor in line mode only; screen 
mode is not compatible with the Tektronix 4014. 

HELP 

Presents a list of subjects and command names for 
which information is available. Provides descrip- 
tions of commands, command formats, and command 
qual i f i ers. 

LOGOUT 

Logs the current user out of the system. 

PURGE filename. typ 

Deletes all but the highest version number of 
file f 1 lename . typ. 


*VAX is a trademark of Digital Equipment Corporation. 
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Table *4.1 Rudimentary Commands (Continued) 


Command 

Description 

SET PR0T-(0:RWED) f 1 lename. typ;n 
SET PROT-(OrRWE) fi lename. typ;n 

and 

change the protection codes for a 
specific file. In the first case 
delete access is allowed; the second 
form protects a file from accidental 
deletion. To determine the current 
protection of all files onter 
DIR/PROT. 


*♦.1.2 SOME MANUALS 
VAX/VMS Primer 

VAX-11 Information Directory and Index 
VAX/VMS Cornmond Language User's Guide 
EDT Editor Manual 

VAX-11 FORTRAN Language Reference Manual 
VAX-11 FORTRAN User's Guide 


*♦.2 DESCRIPTION OF FILES 


Table *4,2. File Descriptions 


Fi le 

Descript i on 

ALLOYS.DAT 

Alloy data base used by models 1 and 2 for 
binary-systems. 

ALL0YS3.DAT 

Alloy data base used by model 3. 

BELL. FOR 

Subroutine which rings terminal bell without 

BELL. OBJ 

moving cursor. 

FTN.COM 

Command file for compllinq. 

GRAPHICS. OLB 

Tektronix software object library. 

ICARD.DAT 

Input card image read by all programs. 

LNK.COM 

Comma i d file for linking. 

Ml .CMN 

COMMON blocks for Ml. 

Ml. EXE 

Executable image for Ml. 

Ml .FOR 

FORTRAN source for Ml. 

Ml M. CMN 

COMMON blocks for HIM. 

Ml M. EXE 

Executable image for M1M. 

Ml M. FOR 

FORTRAN source for M1M. 


* 4-2 







Table 4.2. File Descriptions (Continued) 


FI le 

Description 

M2.CMN 

COMMON blocks for M2. 

M2. EXE 

Executable image for M2. 

M2. FOR 

FORTRAN source for M2. 

M2M.CMN 

COMMON blocks for M2M. 

M2M.EXE 

Executable image for M2M. 

M2M.F0R 

FORTRAN source for M2M. 

M3.CMN 

COMMON blocks for M3. 

M3. EXE 

Executable image for M3. 

M3. FOR 

FORTRAN source for M3. 

MALLOYS. DAT 

Alloy data base for multicomponent-alloy systems. 

MPS.COM 

Command File for executing models 1 and 2. 

MPS3.COM 

Command file for executing model 3. 

PLOTIO.OLB 

Tektronix software object library. 

SPARSE. CMN 

Sparse matrix routines used by moaeis. 

SPARSE. FOR 


SPARSE. OLB 



4.3 HOW 70 RUN OR MODIFY THE PROGRAMS 

To run either the steady-staio model or the unsteady-state model, enter the 
command: 


BMPS name where "name" is one of 

1 binary steady-state model, 

1M multicomponent steady-state model, 

2 binary unsteady-state model, or 

2M multicomponent unsteady-state model. 

Specific operating instructions for each model are in the July 1983 83HV004) 

report. 

To run the coupled-energy model, enter the command: 
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«MPS3 





Specific operating instructions for this model are in the January 198*4 
(8*4HV001) report. 

Two command files, FTh.COM and INK. Cun, have been written to simplify 
modi'ication of the programs. This setup Is intended only for infrequent, 
minor program nodi f icat i ons . 

1. Use s sys.w.n text editor to make your changes in 
Ml. FOR, HIM. FOR, M2. FOR, or M2M. FOR. Note that 
the common blocks for each program are sto r ed in 
the corresponding .CMN file. 

2. Compile the program by using the command 

?FTN program options 

where 'program' is the program name, Ml, M1M, M2, 
cr M2M, and 'options' is any valid DCL FORTRAN 
qua'ifier such as /0°T I H I ZC . 'Options' must begin 
with a /. You can emit 'options' if you want to 
use the FORTRAN default values. The original setup 
was compiled with the default values. 

3. Link the program by using the command 

®LNK program options 

where 'program' is the pregram name as in step 2, 
and 'options' is any valid DCL LINK qualifier such 
as /DEBUG. 'Options' must begin with a /. Yoj can 
omit 'options' if you want to use *.he '.INK default 
values. The original setup was linked with the 
default values. 

After linking, it is recommended that you conserve 
file space by DELEting the file 'program' . OBJ arid 
by PURGing the old version of ' orr^ram ' . EXE. You may 
need to change the file protection codes in order to 
do this. See Section 1.2. 

* 4 . *4 MAGNET I Z TAPES 

These instructions tell how to use tne VAX/VMS system utility called BACKUP 
to write or read a taoe that contains your entire directory. It is the 
simplest way to backup your files or to transport them to another VAX, but o 
tape written by BACKUP cannot be read by any other type of computer. 


Before you can write or read a magnetic tape, you nust have one mounted on a 
tape drive. Get in touch with the computer operations people, and orovide 
them with a tape. The term "MAG TAPE" in the commands below is a logical 


equivalence name for the tape drive. The computer operations people can tell 
you the actual name to use. 

).4.1 HOW TO WRITE A BACKUP TA D E FOR BACKUP OR TRANSPORT 

allocate mag_fape 

MOUNT/FOREIGN MAG_TAPE 

BACKUP t JOHNSTON] MAG_TAPE : SEND. BCK/VER I FY/DENS I TY«l600/REWI ND 
DISMOUNT MAG_TAPE 
DEALLOCATE MAG_TAPE 

4.4 2 HOW TO READ A BACKUP TAPE 

ALLOCATE MAG_TAPE 
MOUNT/FOREIGN MAG_TAPE 
BACKUP MAG_TAPE [ ] 

Dismount mag_tape 

DEALLOCATE MAG TAPE 
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APPENDIX A 


DATA FOR MAR-M246(Hf) 


A. 1 NOMINAL COMPOSITION 


Element wt. pet. Element wt. pet. 


Cr 

9.0 


Mo 

2.5 


Co 

10.0 


Ta 

1.5 


A1 

5.5 


C 

0.15 


Ti 

1.5 


Hf 

1.5 


W 

10.0 


Ni 

bal. 


PARTITION 

RATIOS 





Element 

Ratio, k 



Source 

Cr 


0.9 


Se 1 lamuthu (6) 

in MAR-M200; Wang 





NASAI R-100 


Co 


1.1 


Jeanf i Is et al . 

(8) in Waspal loy 

A1 


1.1 


Sel lamuthu (6) 

in MAR-M200 ( Hf ) 

Ti 


0.95 


Sel lamuthu (6) 

in MAR-M200(Hf) 

W 


1.2 


Sel lamuthu (6) 

in MAR-M2C0(Hf) 

Mo 


0.65 


Wang (7) in NASAI R-100 

Ta 


0.97 


El 1 iott et al . 

(9) in Ni-Ta bi na 

C 


C.25 


Jeanf i Is et al . 

(3) in Waspa 1 loy 

Hf 


0.23 


Sel lamuthu (6) 

in MAR-M200 ( Hf ) 


A. 3 TEMPERATURE DURING SOLIDIFICATION 

The temperature during solidification is assumed to vary linearly 


compos i t ion 
T - T 

o 


according to 
n 

* I ( 3T/ 3C . ) (C 

M 


- C .) 
01 


(A. 1 ) 


where 


n 

T 


T 

0 
C. 

1 

c . 

oi 

( 3T/3C . ) 


number of solute elements, 

temperature (liquidus temperature of the interdendri t ic 
liquid), °C, 

liquidus temperature of the reference alloy, °C, 
composition of solute i, wt. pet., 

composition of solute i in the reference alloy, wt. pet., and 
change in the liquidus temperature due to solute i, °C/wt. pet 


A- 1 


wi th 


As in the previous report (3), values of (3T/9C.) for most of the solutes were 
obtained from the binary phase diagrams (Ni - X) in Elliot et al. (9). For 
hafnium, (3T/3C^) was deduced from Sellamuthu's data on MAR-M200 containing 
0-2.66 wt. pet. Hf (6). 


Element 

C . 

Ol 

wt. pet. 

oT/3C. 

°C/wt. pet. 

Cr 

9.0 

-2.2 

Co 

10.0 

0.0 

Al 

5.5 

-3.2 

Ti 

1.5 

-7.7 

W 

10.0 

+3.5 

Mo 

22.5 

-1.0 

Ta 

1.5 

-61.0 

C 

0.15 

-2.5 

Hf 

1.5 

-10. A 


For MAR-M2A6, the liquidus temperature (T ) is 1360°C (10); therefore, 

T - 1360 -2.2 (C Cr - 9) - 3.2 (C A) - 5.5) 

-7.7 (C T . - 1.5) ♦ 3.5 (C w - 10) 

-2.5 (C Ta - 1.5) - 1.0 (C Mq - 2.5) 

- 61 (c c - 0.15) - 10 . A C Hf (A. 2) 

To apply Eq. (A. 2) during solidification, the composition of each element in 
the interdentritic liquid is computed and Eq. (A. 2) is assumed to apply until 
the temperature drops to 1230°C. At 1230°C, the remaining i nterdendr i t i c 
liquid solidifies as eutectic solid. 

The temperature of 1230°C was selected after consulting Sellamuthu (6) who 
reports that MAR-M200 with 1.5 wt. pet. Hf solidifies with 0.10 fraction 
eutectic (in the absence of macrosegregation). Then by assuming complete 
diffusion of carbon in the solid and no diffusion of the substitutional 
elements, a Schei 1-type calculation shows that the fraction of i nterdendr i t i c 
liquid is 0.10 at a temperature of 1230°C. 


0 

l 
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A. 4 DENSITY OF THE LIQUID 

From the previous report (3), the change in density attributed to each of the 
solutes (excluding Hf) is given. The effect of Hf can be deduced by making 
use of Table 21 in Sellamutha (6); however, his data For 0.86 wt. pet. Hf 
appear to be in error and are not used. 

For MAR-M200 (without Hf), the density is 7-56 g • cm~^ at 1392°C. Using the 
melt compositions as the reference compositions, given by Sellamuthu (6), we 
have 

P L - 7.56 - 0.013 (C C( . - 8.08) + 0.001 (C Cq - 9-31) 

- 0.11 (C A] - 4. 85) - 0.035 (C T . - 1.88) 

+ 0.033- + 0.042 (C u - .2.76) 

Mo w 

- 0.72 (C c - 0.15) + 0.012 (C Nb - 1.05^ 

+ 0.035 C Ja + 0.052 C Hf 

- 0.00125 (T - 1392) (A. 3) 


Here is the density of the interdendritic liquid (g • cm ^) and T is the 
temperature calculated with Eq. (A. 2). The coefficients of the composition 
terms are the respective values of (3p /3 C . ) and the coefficient of the 
temperature term is (3p / 3T) . 


A. 5 DENSITY OF THE SOLID 

The density of the solid during solidification was computed by using the 

scheme presented in Appendix A of the previous report (3). In this report, 

however, some of the coefficients In the equation for the lattice parameter 

have been changed. As before, 
n 

a - a^. + £ (da/dX.) X. (A. 4) 


where 

o 

a N . - 3.5236 A 

a * lattice parameter (cube-edge length in f.c.c. Ni), and 
Xj - atom fraction of solute i. 
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The coefficients (da/dX.) are from the following sources: 


o 


Element da/dX., A 


Cr 

0.105 

Co 

0.019 

A! 

0.186 

Ti 

0.337 

Mo 

0.435 

1 

0.412 

C 

0.065 

Nb 

0.645 

Ta 

1.04 

Hf 

0.75 


The value for hafnium is assumed to be 
two elements are chemically similar an< 


Source 

Loom is ( 1 1 ) 

Loomis (11) 

Loom is (11) 

Loomis (11) 

Loomis (11) 

Loomis (11) 

Structure Reports (12) 

Loomis (11) 

Calc, from at. radi i 
Kornilov and Snetkov (13) 

equal to that for zirconium since these 
have equal atomic radii. 


With the lattice parameter calculated according to Eq. (A. 4), the densities of 
the nickel-rich phase (y) at 20°C were calculated. Compositions selected were 
those corresponding to Scheil-type solidification. 


The densities of y at the elevated temperatures in the solidification range 
were then calculated by using the expansion coefficients presented in the 
previous report (3)* 


The density of the eutectic-solid is not available so a calculation using a 
number of simplifying assumptions was made. The first ass'imption is that 
eutectic-solidification occurs isothermal ly at 1230°C. From the work of 
Sellamuthu (6), the eutectic comprises Y+Y 1 and carbides in MAR-M200 with 2 
wt. pet. Hf. To estimate the portion of the eutectic which is carbide, we 

assume that the carbon rsr ning in the interdendr i t ic liquid reacts 
completely to form HfC at 1230°C. According to this simple calculation, the 
eutectic comprises approximately 2.5 pet. HfC and the balance is a mixture of 
Y/y'. Using the composition of the Y/y 1 , the calculated density is 8.052 
g/ern^ at I230°C. The small amount of carbide in the eutectic is ignored; 
therefore, the density of the eutectic solid is 8.052 g • cm 



A. 6 SUMMARY OF DENSITIES 


Figure A. 1 summarizes the calculated densities of the liquid- and solid-phases 
during solidification. Figure A. 1 gives the densities versus temperature; it 
shows that the behavior of MAR-M2**6 (Hf) is substantially different than 
MAR-M246 without Hf. With 1.5 wt. pet. Hf, the total density change of the 
i nterdendr i t Ic liquid increases by only pet., whereas with no Hf the 
density of the interdendritic liquid decreases by 16 pet. (see Fig. A. 2 in the 
previous report (3)K With hafnium, the density of the solid is approximately 
the same as that for the alloy without hafnium. 

The effect of hafnium on the density is attributed to its direct effect on 
increasing density and, also, to its effect on the partition ratios of some of 
the other elements (especially A1 and Ti). 

For the calculation of convection and macrosegregation, it is recommended that 
Eqs. (A. 2) and (A. 3) be used to calculate the density of the liquid (pi_). The 
recommended values for the solid-densities are P s 3 7.82 g • cm"3 and 
P SE * 8.05 g * cm"^. 
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Figure A. 1 . 


Densities of Liquid- and 
During Solidification of 


Solid-Phases and Fraction 
MAR-M2^6( Hf ) 


Solid 
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TEMPERATURE. 


APPENDIX B 


DATA FOR PVA 1480 (ALLOY 454) 


B. 1 NOMINAL COMPOSITION 


Element wt. pe t. 


Cr 

10.0 


Co 

5.0 


W 

4.0 


Ti 

1.5 


Ta 

12.0 


Al 

5.0 


PARTITION 

RAT 1 OS 


Element 

Ratio, k 

Source 

Cr 

0.9 

Sellamuthu (6) in MAR-M200 

Co 

1.1 

Jeanfils et al. (8) in Waspalioy 

W 

1.2 

Sellamuthu (6) in MAR-M200 

Ti 

0.55 

Sellamuthu (6) in MAR-M200 

Ta 

0.70 

Kadalbal et al. ( 1 4 ) in Ni-Al-Ta 

Al 

0.95 

Kadalbal et al. (14) in Ni-Al-Ta 


B.3 TEMPERATURE DURING SOLIDIFICATION 

The temperature during solidification is assumed to vary linearly with 
composition and temperature according to Eq. (A.l) and the coefficients of 
Section A. 2 apply. A reference temperature of 1412°C for N i -5 - 8 1 wt. pet. 
A1 - 15.18 wt. pet. Ta is used (14). The resulting equation is: 

T - 1412 - 2.2 C Cr - 3.2 (C A] - 5.81) 

- 7.7 C Tj + 3-5 C y - 2.5 (C Ja - 15.18) (B. 1 ) 

Equation (B.l) gives a liquidus temperature of 1 40 ^ 0 ^ jn d applies to 1290°C 
when eutectic solidification occurs. The latter temperature was selected 
because Gel 1 et al. (15) reported that incipient melting in Alloy 454 occurs 
at 1288°C. 
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B 4 DENSITY OF THE LIQUID 

! t is assumed that Eq. (A. 3) applies; so that for no molybdenum, carbon, 
niobium, and hafnium, the density of the i nterdendr i t ic liquid is 

P L - 7.66 - 0.013 (0 Cr - 8.08) + 0.001 (C Cq - 9.31) (B.2) 

- 0.11 (C A1 - 4.85) - 0.035 (C T| - 1.88) 

+ 0.04?. (C w - 12.76) + 0.035 C Tg - 0.00125 (T - 139?) 

B.5 DENSITY OF THE SOLID 

Equation (A. 4) and the coefficients listed in Section A. 5 are used for PWA 
1480. The densities of Y at the elevated temperatures were calculated by cal- 
culating compositions for Schei 1-type solidification with the partition ratios 
of Section B.2 and then by using the expansion coefficients in the previous 
report (3). 

I 

Because PWA 1 480 contains no carbide, the eutectic comprises Y -Y and a 
small amount of Y-Y 1 - <5 which solidifies over a small temperature ranr*. 
(14). As discussed in Section B.4, the solidification of the eutectic is pre- 
sumed to be isothermal and at 1290°C. 

B.6 SUMMARY OF DENSITIES 

Figure B.1 summarizes the calculated densities of the liquid- and solid-phases 
during solidif cation as well as the temperature versus fraction solid. 

For convection and macrosegregation calculations, Eqs. (B.1) and (B.2) should 
be used to calculate the density of the liquid during solidification. For the 
primary solid-phcse, ° s - 7.77 q • cm ^ and P^ - 7.86 g • cm 
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Figure 6.1. Densities of Liquid- and Solid-Phases and Fraction Solid 
During Solidification of PWA 1^80 (All ov 
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APPENDIX C 


FORMATION OF POROSITY IN NICKEL-BASE ALLOYS 

Ouring solidification, hydrogen segregates to the i nterd ^ndr I t i c liquid and, 
if porosity forms, reacts according to 

iH 2 (g) - H. (C.1) 

in which H Indicates that hydrogen is dissolved. 

Similar reactions can he written for nitrogen and for oxygen, but the partial 
pressures of these two elements are negligible when compared to the partial 
pressure of hydrogen. In the case of nitrogen, almost all of the alloy 
elements in nickel-base superalloys have the effect of strongly reducing the 
activity of nitrogen (so that the partial pressure of nitrogen is negligible) 
and/or the effect of forming nitride compounds. Similarly, the reactive 
elements in nickel-base superalloys combine with the oxygen to form oxides. 

In steel, it is well known that carbon can react with oxygen to form carbon 
monoxide during solidification: 

C + 0 - C0(g) . (C.2) 

However, in the presence of strong "deox i d i zers" reaction (C.2) does not occur 
in steel. The reactive elements in nickel-base superalloys-, in this regard, 
are strong "deoxidizers" so that the evolution of carbon monoxide during 
solidification does not occur. 

The free energy of reaction (C.1) in liquid nickel is ( 1 6 ) 

AG° - *«800 + 8.36T ( C. ? ) 

in which the standard state for dissolved hydrogen is taken to t _ 1 wt. pet. 
Because 

AG° - -RT in a H /p H i (C.A) 

2 

then 

P H - a* exp (3.A15 ♦ '»831/T) (C.5) 


C-l 


where 


p u ■ partial pressure of hydrogen, atm. 

H 2 

T - temperature, K 

-1 

R ■ gas constant, 1.987 cal(gm - atom) K 

and 

a u « activity of hydrogen dissolved in the liquid nickel. 


The activity of hydrogen is 
5 H " f H C H 


(C.6) 


in which f^ is the activity coefficient of hydrogen and is the concentra- 
tion of hydrogen in wt. pet. In multicomponent superalloys, the activity 
coefficient is not known so we resort to a technique, which is strictly 
applicable to dilute alloys, to estimate f^. The alloy elements affect ? H 

according to 

n i 

log f H “ l c : ( c - 7' 

J J 

in which er* are the interaction coefficients of each n alloy element. 

H 


Sigworth et al. 
given below. 

Element 


(16) have summarized many values of e^. some of 

M 

Element 


wnicn are 


'H 


e H 


Al 

0.014 

Ti 

0.013 

Co 

0.0031 

Ta 

0.011 

Cr 

0.0036 

Hf 

0.011 

Mo 

0.011 

C 

G.C 65 

W 

0.011 




Sigworth et al. (lb' presented no data for Ti, Ta, and Hf, so estimates were 
made based upon the proximity of thece elements to others, for which is 
known, in the periodic table. The value for carbon (i.e., e^) is assumed to 

n 

be that for e^ in liquid iron (17). 

H 


The partitioning of Sy:'rogen between the liquid- and solid-phases in nickel- 
base superalloys is not known, but some data are available for nickel, nickel- 
cobalt alloys, and nickel-iron alleys (1^). To extrapolate these data to 
obtain values of partition ratios, solubility data were plotted as 


C-2 



in s vs. 1/T as indicated in Fig. C.1 where s is the solubility of hydrogen in 
equilibrium with 1 atm. of ^(g). A linear relationship is expected provided 
AG° is of the form given in Eq. (C.3) and the activity coefficient is 
independent of temperature. 

From Fig. C.1, the partition ratio for nydrogen was determined at 1200OC, 
13U0°C, and 1400°C for nickel and N i -20 wt. pet. Fe and at 1400°C for Ni-20 
wt. pet. Co. Values ranged from 0.42 to 0.46. A value of k u = 0.44 is recom- 
mended. 

Sample calculations were done for hydrogen dissolved in PWA 1480 with the 
properties given in Appendix B. The substitutional alloy elements (Cr, Co, W, 
Ti, Ta, and A1 ) were assumed to partition according to Schei 1-type solidifica- 
tion and dissolved hydrogen was assumed to partition with uniform concentra- 
tions within both phases. 

The scheme employed for the calculations follows: 

1. Calculate the composition of the interdendr i t ic liquid; 

2. Calculate the activity of dissolved hydrogen in the i nterdendr i t i c 
liquid with Eqs (C.7) and (C.6); and 

3. Calculate the partial pressure of ^(g) with Eq. ( C . 5 ) • 

Results of such calculations are shown in Fig. C.2. With a melt containing 5 
ppm of dissolved hydrogen, the pressure generated is less than 1 atm. until 
near the erd of solidification whereas with 8 and 10 ppm, the pressure gen- 
erated is substantially greater. Presumably, increasing amounts of micro- 
porosity would be expec’.ed in, the order of, nelts containing 5, 8, and 10 
ppm. 

To nucleate a gas bubble within the interdendr i t ic liquid, the internal hydro- 
gen pressure must overcome the local pressure within the solid-liquid region 
and the added pressi re due to surface tension. 
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Spec i f ica I ly , 


p u > p + 2 Y/R fC.8) 

. H 2 " 
where 

p * pressure inside the bubble, 

H 2 

p - local pressure, 

Y = surface tension, and 

R = radius of gas bubble. 

For nickel, y ■ 1778 mN • m ^ (19). As an example, assume that p = 1 atm 

(1.013 N • m~^) and p,, * 1.2 atm (1.216 N • m ^). From Ea. (C.8), it follows 

ru 

2 

that a bubble is stable provided R L 1.65 x 10 m or 1 65 ym. A typical pri- 
mary dendrite arm spacing in directionally solidified superalloys is 300 ym; 
so, presumably this calculation indicates that microporosity would form. 

The use of Eq. (C.8) gives an overestimate of p u required to produce a stable 

H 2 

bubble because there is no account of the presence of the solid-liquid inter- 
face that would provide an active site for bubble nucleation. 



NOTE: Solubility (S) has units of cnr of H ( g ) at STP/100 g of metal 


Figure C. 1 . 


Solubility of Hydrogen in N,„*el and N i -20 wt. pet. Co and 
N i —20 wt. pet. Fe Alloys 


r-t; 






PARTIAL PRESSURE OF H~ ATM. 



NOTE: Numbers on the curves refer to the concentration of dissolved hydrogen 
in the melt. 


Figure C.2. Partial Pressure of H„(g) Generated in Interdendr i t ic 
Liquid During Solidification of PWA 1480 





APPENDIX D 


SURFACE TENSION IN NICKEL-BASE ALLOYS 

A relation derived by Guggenheim (20) for ideal binary solutions was applied 
by Speiser and Spretnak (21) to estimate the surface excess at a grain bound- 
ary; the relat ! on is 

exp( -Y s/kT) = N 1 exp(-Y 1 s/KT) + N 2 exp( -Y 2 s/kT) (D.l) 

where 

Y = surface tension, 

s = surface area per atom, 

T = absolute temperature, 

< = Boltzmann's constant, and 

N^, N 2 = atom fractions. 

Speiser (22) has suggested that, lacking data, Eq. (D.l) could be applied to a 
multicomponent alloy, and it would be better to account for the individual 

atomic surface areas of each component. Then Eq. (D.l) becomes 

v 

exp(-Ys/*T) = . N. exp( -y. s . /<T) (D.2) 

where s. is the surface area per atom, and s is the weighted average (based 
upon the atom fractions of all components). 

The surface tension of each component is given by 

y = Y q + (T-T q ) (dY/dT) (D.3) 

where Y q is the surface tension of the component at its melting point T . 

Equations (D.2) and (D.3), together, can be used to calculate the surface 
tension of the interdendr i t ic liquid during solidification; the calculation 
accounts for the change in composition and temperature of the i nterdendr i t i c 
liquid. Using the data and relations in Appendix A and Appendix B, the compo- 
sition and temperature of the i nterdendr i t i c liquid in MAR-M2A6(Hf) and PWA 
1480 were computed assuming Scheil-type solidification. Then Eqs. (0.2) and 
(D.3) were used to calculate the surface tensions shown as the two lower curves 
in Fig. D.l. 
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The atomic surface areas were computed from the atomic radii given in Brandes 
(23); surface tensions and values of (dY/dT) were also taken from Brandes 
(19). These data are summarized below for the elements in MAR-M2*46(Hf ) and 
PWA 1*480. 


Element 

T 0 , °C 

Y q , dyn'em -1 

dY/dT, 

-1 r -1 
dyn’em *C 

S i ’ 

cm 2 x10 16 

Cr 

1875 

1700 

-0.32 

5. *48 

Co 

1*»93 

1873 

-0.*49 

3.95 

A1 

660 

91*4 

-0.35 

6.-42 

Ti 

1685 

lA«;o 

-0.26 

6.80 

W 

3377 

3 

-0.29 

6.25 

Mo 

2607 

22 3 0 

-0.30 

6.15 

Ta 

2977 

2150 

-0.25 

6.80 

C 

-- 


— 

1.86 

Hf 

19**3 

1630 

-0.21 

7.95 

Ni 

1*45** 

1778 

-0.38 

*4.90 


Carbon does not affect the surface tension of iron (2*0; it is, likewise, 
assumed not to affect the surface tension of these alloys except by the small 
effect when the average value of s for the alloy is computed for use in Eq. 
(D.2) . 

Figure D. 1 shows that the surface tensions of the interdendr i t ic liquid in 
MAR-M2*46(Hf ) and PWA 1 480 are approximately 5-10 pet. less than the surface 
tension for pure nickel, curves a and b. Curves a and b show a slight 
variation in surface tension only because the temperature during 
solidification varies whereas the curves for PWA 1*480 and MAR-M2*:6( Hf ) reflect 
both variations: temperature and composition. Clearly, the effect of 
composition predominates. 

Because the surface tensions of the alloys vary only slightly during solidifi- 
cation, it is recommended that constant values be used to predict the forma- 
tion of porosity during solidification. For MAR-M2*46( Hf ) , the recommended 
value for the surface tension is 1700 dyn*cm ^ . For PWA 1*480 , it is 1710 
dyn*cm _1 . 
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FRACTION SOLID 


Curves a and b show the effect of temperature alone on the surface 
tension of nickel. 


Figure D. 1 . 


Surface Tension of I nterdendr i t i c Liquid huring 
Solidification of T wo Nickel-Base Superalloys 






APPENDIX E 


FORMATION OF POROSITY IN ALUMINUM-COPPER ALLOYS 

The hydrogen reaction during pore formation in aluminum-copper alloys is the 
same as for pore formation in nickel -base alloys: 

iH 2 (g) - H , (E. 1 ) 

in which ^ indicates that hydrogen is dissolved. The expression for the 
partial pressure of hydrogen is then 
2 

P H 2 = a H exp( 5.96A + 12508. /T ) (E. 2) 

where 

p H ■ partial pressure of hydrogen, atm., 

2 

T = temperature, K, 

and a H = activity of hydrogen dissolved in the liquid aluminum. 

The activity of hydrogen is 

a H ’ f H C H < E '3) 

in which f^ is the activity coefficient of hydrogen and C H is the concen- 
tration of hydrogen in weight percent. The effect of copper on the activity 

coefficient of hydrogen in liquid aluminum is 

In f H = 0.07A8 C L - 0.001A2 C* at 750°C, and (E.A) 

In f H - 0.053^ C L - 0.0008A8 C* at 950°C, (E.5) 

where is the weight percent copper in the fluid. To obtain the activity 
coefficient at temperatures in the solidification range, 5^8°C to 660°C, Eq. 
( E . 4 ) and (E.5) are extrapolated linearly in the reciprocal of temperature: 



(E.6) 


f 4 - A C L + B c[ at T °C, 

where A - + (A 2 -A ] )(l/T-1/T 1 )/( 1/T 2 -1/T 1 ), (E.7) 

B - B, + (B 2 -B.)(1/T-1/T 1 )/(1/T 2 -1/T 1 ), (E.8) 

and the subscripts 1 and 2 reference values at 750°C and 950°C, respectively. 

The concentration of hydrogen in the interdendr i t ic liquid is given by the 
lever rule, 

C H = C HI /(g l + k H g S ) » (E ’ 9) 

where the partition ratio for hydrogen, k^, is 0.0555 and C H| is the initial 
weight percent of dissolved hydrogen in the molten metal. 
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APPENDIX F 


SURFACE TENSION IN ALUMINUM-COPPER ALLOYS 


The surface tension of the interdendrl t Ic liquid during solidification of 
aluminum-copper ail >ys can be calculated In the same manner as for nickel-base 
alloys (Appendix D). The calculation accounts for the effects of variations 
in composition and temperature of the interdendr i tic liquid on the surface 
tension. Surface tension, temperature, and liquid composition of Al-4.5 
wt.pct.Cu are shown in Table F. 1 as functions of the volume fraction solid for 
a Scheil solidification calculation. 


Table F.l. Scheil Solidification Calculation Results 


Vo 1 ume 
Fract ion 
Solid 



Liquid 

Temperature 

Surface 

Compos i t ion 

(°c) 

Tension 

(wt.pct. ) 


(dyn/cm) 

4.5 

645 

949 

5.65 

641 

958 

7.07 

636 

969 

9.62 

628 

987 

15.78 

607 

1029 

30.34 

558 

1 118 


Because the surface tension varies significantly during solidification, it 
should be written as a function of temperature or composition. In terms of 
temperature, 


o - 85**.** + 2.555 T - 0.003734 T 2 


(F.l) 


where o is in dyn/cm and T is in °C. In terms of composition, 


o - 915.0 + 7.828 C L - 0.03736 C* 


(F. 2) 


with o in dyn/cm and C L in weight percent Cu in the interdendr i t i c liquid. 
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